


I . 

mm WWW 1 


Prepared by: Win 
Ittcal Mechanics 


D ANALYSE 


$ AND SPACE AD 


NED SPACECRAFT CENTER 


mUSim, TEXAS 


June IS - 


ImOSSSON NUMBER) 


“ 37 (PAGis) ^ „ 

/ /v? y A //& 

X lZIA; 

{NASA CR OrTmX OR AD NUMBtR) 


(CODE) 


(CATEGORY) 



MSC INTERNAL NOTE NO. 66-FM-55 


PROJECT APOLLC 

THE GENERALIZED FORWARD ITERATOR 

Prepared by: William E. Moore 
Analytical Mechanics Associates, Inc. 


June 15, 1966 


MISSION PLANNING AND ANALYSIS DIVISION 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
MANNED SPACECRAFT CENTER 
HOUSTON, TEXAS 


Approved: 



Mission Analysis Branch 


Approved: 


John#. Mayer, thief * 

Mission Planning and Analysis Division 



1 


THE GENERALIZED FORWARD ITERATOR 
by: William E. Moore 


SUMMARY 

A computer program has been written which constructs optimized 
Apollo missions. It is highly versatile; it will produce a scan of 
several related missions or an accurate trajectory of reference qua- 
lity. Speed and accuracy are purely dependent on the mathematical 
model chosen for each particular case. This model could consist of 
an array of two-body approximations, or of a precise numerical inte- 
gration of the equations of motion going by the Eneke method when 
applicable. The program is extremely flexible with regard to the 
mission that can be constructed, and the extent to which it can be 
optimized. This note indicates how the input describes the mission 
and the method of computation to the computer. It also presents the 
various mathematical trajectory models that are available. 


INTRODUCTION 

This Apollo mission design program has been designed to furnish 
the analyst with a reliable tool for the design of a great variety 
of missions. It is fast enough to be capable of producing a "scan” 

when no great accuracy is required, as in preliminary mission design, 
and accurate enough to produce a trajectory of reference quality 
when required. It is general and flexible enough to permit the ana- 
lyst to assess the effects of many different mission configurations 
and the introduction or omission of numerous constraints. 

The major constituents of the program are (l) input and initial- 
zation, (2) first guess generation, (3) trajectory, mission parameter 
and constraint computation, (4) iterative parameter correction for 
the satisfaction of constraints and eventual mission optimization 
procedure, and (5) program output. These constituent parts are 
described in detail in the following sections. 
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PROGRAM INPUT AND INITIALIZATION 

The operation of the program is completely controlled by a set 
of input data and switches. Each input number has associated with 
it on a card its location in the full input array. This location 
appears on the card and the input processor then stores this number 
in its proper lc nation. Thus, it is necessary only to load those 
input numbers which are needed and which differ from previously loaded 
value. Physical quantities may be read in a number of different 
units. The unit is indicated by a code letter on the input card, 
which enables the processor to apply the proper scale factor for 
conversion to internal units. The input processor further has the 
capability of using a previously loaded number, of incrementing it by 
a specified amount, or of using a specified multiple or fraction of 
the present value. 

A first guess for the independent variables obtained previously 
may be retained, be incremented or else be recomputed for every case. 
In addition, the converged answer of the previous case may be used 
as a first guess. 

The input array further contains a set of switches which substan- 
tially control the flow of the program. 

One set of these describes the situation at the beginning of the 
trajectory, and thus indicates which maneuvers are still to be perfor- 
med. For example, the trajectory might start at the launch site, or 
at a later time, such as same position and velocity in earth parking 
orbit, or along the trana lunar leg, in lunar parking orbit, or along 
the transeartb leg. 

Another set of switches describes the independent variables of 
the search; i.e., the variables to be determined, and the constraints 
to be satisfied. Table II gives a list of all the inputs to the 
program. 

During initialization all these switches are used to set up a 
path through the trajectory computer, so that only the trajectory 
segments needed in the problem at hand are computed, and all neces- 
sary constraints evaluated. It further guarantees that all neces- 
sary partial derivatives are computed. Simultaneously, in addition, 
all these switches are checked for consistancy. If an obvious incon- 
sistency is detected, the program proceeds to the next case. Because 
of the great generality of the program, a complete consistency check 
is not feasible. Therefore, even if the input passes all the 
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• e*>isistency tests, this "consistent" trajectory is not guaranteed to 

represent a reasonable problem, but as much checking of the switches 
as possible is done. 

A third set of switches controls mode of computation of each of 
the various sub-arcs in the trajectory. The alternatives here are 
computation by numerical integration or by conic approx i at ions for 
the coast phases and idealized arcs for the powered portions. 

The choice here will depend on the accuracy required as well as 
on the computer time available. 


FIRST GUESS GENERATION 


In order to insure the most efficient behavior of the iteration 
scheme, in some cases, first guesses for one or more of the indepen- 
dent variables may be computed. No effort has been made to fill up the 
program with all the known direct methods for getting these values; 
rather only those which make the trajectory computer program using or 
the iterator operation more efficient are included. 

One can choose to compute first guesses for time of launch and 
time in earth parking orbit if starting at the launch site. Falling 
out, also, is the velocity increment at translunar injection. These 
three values are based on the desired translunar flight time, or a 
nominal time if a free-return trajectory is desired. From this value 
and the base time, the position of the moon at some point after arri- 
val is computed as a target. Now the time of launch is chosen so that 
the plane of the earth parking orbit contains this target point. Then 
the time in that orbit is adjusted so that the translunar trajectory 
passes through the target point. These times are close enough to the 
correct ones to be quickly adjusted by the parameter search procedure. 

In the case where the parking orbit is already fixed, the proce- 
dure is slightly modified in that the translunar trajectory — still 
assumed to be in the same plane as the parking orbit — commonly will 
be chosen as close as possible to the target point, but will not, in 
general, contain it. In this case, too, the search procedure quickly 
adjusts the first guess to their correct value. 

Another case involves a point after translunar injection, but 
not on the nominal translunar trajectory. In this situation, the 
most efficient procedure turns out to be to apply the search technique 
in a "backward" manner first. This consists of searching for a time 
of pericynthion and a pericynthion velocity vector, given a fixed 
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pericythion position, which will lie on a trajectory containing the 
original given point. At the latter point, the computed velocity 
vector provides first guesses for the maneuver to be executed to return 
to nominal condition. 

Finally, for each of the major burn phases, provision is made to 
compute first guesses for the guidance parameters. To use this provision, 
the analyst must furnish a vector approximating the state at the end 
of burning. 


TRAJECTORY COMPUTER 

The variables to be determined are, of course, initial and 
control variables; these are called independent variables. The con- 
straints to be satirfied by the mission are represented by parameters 
called dependent variables. They are quantities which can be derived 
from states at significant points along the trajectory. Thus, by 
the trajectory computer we mean the program which transforms the 
independent variables into state variables, follows the values of the 
state variables along the mission path, and transforms the state 
variables into the dependent variables. The transformations are usally 
simple computations based on geometrical or physical relations and may 
cn occasion involve a simple iteration. The trajectory proper - the 
propagation of the state variables along the mission path - is more 
complicated. Computing the state variables accurately is a time-con- 
suming process. Thus, the user has a choice between an accurate but 
slow, and a fast but approximate method of computing each arc. These 
arcs are divided into two classes, according as power (that is, thrust- 
ing) is off or on. The arcs without power are called coasting arcs 
and include earth parking orbit, translunar coast, linear parking orbit, 
free-return coast, and transearth coast. Each of these may be integra- 
ted, if accuracy is desired, by the Encke method with the "universal 
incremental anomaly" as independent variable. Where speed is more 
important, as happens when there are many cases to be run(e.g. in a 
scan of initial or final parameters) the geometrical-dynamical equations 
of two-body motion are used. 

The Encke method is, by, now, reasonably familiar. It consists of 
separating the accelerations into two-parts — one of which gives a 
differential equation solvable in closed form (the "two-body" solution), 
and the other of which is small, (the "perturbations" ) so that, even 
though they must be integrated numerically, the integrations can proceed 
in relatively large steps without the introduction of appreciable 
errors. The programming for the two-body solution if available anyway, 
in connection with the approximate trajectory methods; hence, no extra 
space is used to handle it. Moreover, each integrated arc of the 
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trajectory has a certain termination condition, which is the product 
of an iteration in step size. First guesses for stopping conditions 
follow easily from the two-body formulas. 

Both the two-body and perturbation equations are written in terms 
of the anomaly 3 and derivatives with respect to 3. Doing this allows 
the integration step size to adapt Itself to changes in distance from 
the central attracting body. Furthermore, the computation of the time, 
t, from 3 is decidedly easier than computing £ from t. This is because 
Kepler's equation is transcendental in 3, but linear in t. For a 
description of the equations in terms of 3> see the Appendix. 

The equations for the two-body motion are formulated in terms of 
two different anomalies. The "universal incremental anomaly", 3, and 
the "eccentric incremental anomaly", 9. It would be possible to write 
all the equations in terms of 3 but in this way computing time has 
to be sacrificed to elegance; hence, the "universal" approach has been 
abandoned where necessary. The current method is popularly referred 
to as the "Herrick- Beta" method. 


The formulas are divided into two classes. The first class contains 
procedures for determining the value of the universal anomaly through 
which the position and velocity and time must be propagated. Here 3 
is generally described in terms of the eccentric anomaly, 0. The 
second class contains the procedures for propagating the position and 
velocity vectors, and the time. In all that follows X , X , and t are, 
respectively, the position vector, the velocity vector? anS the time at 
a known point on the trajectory. Let n be the mass parameter of the 
attracting body. 


We start with the second class of formulas, that is, the propaga- 
tion formulas. Let 
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( 1 ) 

( 2 ) 

(3) 


1 = 2 - X2 2 


a r 0 p 



(*) 

(5) 

oft 

2k + i)! (i = 0, 1, 2, 3) (6) 
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X = f X + g X 

o o 


• • • • 

X = f X + g X 

o ® o 


(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 

( 12 ) 

(13) 

( 1*0 

(15) 


Formulas (9), (l4), and (15) perform the propagations, hut formulas 
(7) and (8) are also convenient in some applications. 


The equation for p depend on some specified configuration of the 
state such as periapsis, specified distance, and specified flight-path 
angle. For periapsis, compute ro, v 0 ^, d Q , and ^ as in equation (l) - 
Ik) above. 



If 1 >0, set 
a E = tan 



taking all four quadrants 
into considerations 


(19) 
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If 1 <0, set 

a 

E = log s + c 
e 



(19’) 

( 20 ) 


This arrangement breaks down if 1 = 0, but this case is impossible 
to create in the computer. a 


To determine the ^ for a given distance, it is first required 
that the current distance ro is less than, or equal to, the desired 
distance s„ It is assumed that, in the case of elliptic orbits, 
no more than one revolution is traversed frcm input state to output 
state. With these two assumptions, for each direction, incoming and 
outgoing, there can be no more than one place on the orbit at the distance 


r . 


Thus a switch, 



W, is provided, 
if forward 
if backward 


( 20 ) 


Again formulas (l) - (4) provide d and 1, and (17) and (l8) are 
used to get c and e. 0 a 

Now set C = c (21) 

If "g >0, and C < -1, the distance r is impossible, and calculation 
is suspended; otherwise, compute 


V 11 - C o S 1 

(25) 

- c 8 ! 

(2U) 


Now if 1 > 0 

(25) 

(26) 

with the results allocated between quadrants I and II. 
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(Four- quadrant allocation does this automatically because (23) and 
(24) give non-negative results.) 

If TC <0 

. E 0 - log (C o ♦ S Q ) (25') 

and E -= log (C + S) (26) 

Now |0| =| (sgn d o ) E q - W E| (27) 


represents the absolute value of the total 
eccentric anomaly. Tfous 


= W \B\ 



(a8) 


To determine the 0 for a given flight -path angle we require that 
the input state must be atperiapsis. Formulas (l), (2), and (4) arc 
again used to get r Q and 4. Now the eccentricity, e, is given by 


e = 1 - 


Compute 


= ^| e 2 - U 

frk- 

jifT 7” 


{ 


tan" 1 £ 
C 


s in y 
cos y 

if 1 > 0 
a 

if 1 < 0 
a 

if 1 » 0 ■ 

a 


log C + C if 1 
e a 


< 0 


and finally, 


0 



(29) 

(30) 

(31) 


(32) 


(33) 
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Note that Jn the case of an elliptic orbit, formula (31) allows 
only that occurrence of the flight-path angle which is close to per- 
1 apsis. 

If any arc has to pass from the vicinity of one attracting body 
to the vicinity of another an iterative procedure is used to make 
sure that the arc is a continuous piecing together of two-body arcs. 

At the point where the arcs .loin, both the positions and velocities 
must agree. Hiis "patch" point is defined as the intersection of the 
initial two-body arc with the surface around the moon where the ratio 
of the distance from the moon to the distance from the earth has a 
given value. 

To start the procedure, the initial state is used, by means of 
the distance formula previously given, to find two values of the univer- 
sal anomaly 0 for which the ratio described above will bracket the 
given value. If this is impossible, the precedure ie terminated because 
the initial state is not such as to determine an arc intersecting the 
surface. From these first two values of 0, a third value is obtained 
by linearly interpolating to get the desired ratio. Then these three 
values are used to start a quadratic iteration process, at each stage 

of which the outer two of three values of 0 used are chosen to bracket 
the value of the ratio, and the inner value is the last value computed. 

The process is stopped when this last value changes by a negligible 
amount . 

For the burning arcs in the integrated mode, the thrust accelera- 
tion is incorporated among the terms to be integrated. It is computed 
from the guidance logic appropriate to each maneuver. In the approxi- 
mate mode a geometrical representation of the effects of burning is 
used. That is, the changes in time, distance, velocity, and flign J path 
angle are applied, and the effects of rotation in the previous pVoe 
of motion, and out of it, are incorporated. The changes and amounts 
of rotation have to be calculated beforehand by studying the results 
of the integrations of the burn arcs. Let x 0 , x Q be the position 
and velocity vectors at the beginning of the burn arc. Let Ar, Av, 
and AY be the changes in distance, velocity, and flight-path angle, 
respectively. Let As, AA, be the amounts by which the position vector 
is rotated in the plane of motion, and by which the velocity vector 
is rotated out of the plane of motion, respectively. 

2 

Let r = x .x 
o o 



d 
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b ■ 1*0* x o' 


Wien x, = x cos As + r x - dx sin Ax 
l o o 


• ^ 

x, = x_ cos (As - A v ) + ^O v x o sin (As 

1 o — 

h 


- Ay) 


Let d^ = x^ • x 1 


= (x x . x x f 


Then x^ = x. 


2d i : 


AA 


x„ = — -g- sin ^ + x^ cos AA - X 1 * X 1 sin AA 


- (*„ . = r, and v = (x„ . x„f 


2 2 


IlnaUy '’ I, * ^1 

*3 = X 2 U + ^ 


2 2 ' 


*3 


= (1+^) 


are the position and velocity vector at the end of the burning arc. 

( \x - x |\ 

1 3 o 1 where I is the specific 

^ / 

impulse and g is the acceleration due to gravity at the earth's surface. 


A narrative description of the actual method of computing the 
trajectory follows. The independent and dependent variables are referred 
to as xi and y^, respectively. For example, X25 is the independent 
variable time of launch, and y3 is the dependent variable time of launch. 
Many quantities required by the trajectory computer may be computed, once 
for all trajectories, as soon as the input is available. Thus the reader 
need not assume that every computation shown is performed every time the 
trajectory computer is passed through. Instead, the trajectory computer 
is intended to take care of the effects of changing values of any of the 

Xi . 
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We start at the launch site. Knowing the latitude and longitude 
of launch, and the time of launch, x.-, we know the position of the 
spacecraft at launch. An art if icial velocity vector is made up out 
of the launch azimuth and the circular velocity at the height defined 
by the Ar for ascent. Now an application of the approximate burn 
formulas furnishes the state at the end of ascent (or at the beginning 
of earth parking orbit). In the process the deviations in flight -path 
angle and velocity from circular conditions, X22>and X23, are used. 
Note that there isn't any provision for integrating the ascent trajec- 
tory. 


Hie time and azimuth at launch, y, and y. , are simply transferred 
over as dependent variables. Two-body formulas are now applied to the 
state at the beginning of earth parking orbit to get y^ and y<, the 
perigee and apogee, heights, respectively, of the orbit. If the 
orbit is circular, the initial height is used. Next the state at the 
beginning of earth parking orbit is propagated to the state at the end 
of earth parking orbit, by using the time in earth parking orbit, x_, 

Hi is computation may be either accurate or approximate. 

For the approximate translunar injection change in velocity and 
plane change, x„ and x, Q , are used in the burn formulas to get the 
state at the enaof translunar injection (or the state at the beginning 
of translunar coast). For the accurate translunar injection, the ob- 
solete MIT guidance is used, pending an exact definition of something 
better. In either case, the mass after translunar injection, y_, is 
available . 

The state at the beginning of translunar injection is either 
propagated through the patching iteration, or integrated to pericynthion, 
where the height, yg, and translunar flight time, y^, are available. 

Hie flight-path angle at the beginning of lunar orbit insertion, 
x^4, is used, together with the state at pericynthion, to get the state 
at the beginning of the lunar orbit insertion maneuver. Hiis is done 
either by the approximate formula given above or by integrating accura- 
tely to a state which has the given flight -path angle. At the start 
of lunar orbit int-rtion, the height xg, is available. The coordinates 
of the state are then transformed so that they are referenced to the 
earth-moon plane, and the inclination, latitude, longitude, and azimuth, 
yiO, yn , yi2> and y^x, are computed in that reference frame. 

If free-return constraints pre to be computed, the state at the 
beginning of lunar orbit insertion is propagated back to the earth 
to get a state of perigee, from which the height and inclination, y ^5, 
and y-^g, can be obtained. 
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Next the state at the start of lunar parking orbit Is used to 
compute the effects of the lunar orbit Insertion maneuver. In the 
process, either X13 or xqq and xqg are control parameters for this 
manever. As a result, we get the state at the beginning of lunar 
parking orbit and the height yq8 and mass yqy at that point. If the 
parking orbit is to be non-circular, X9 and x^q, the flight-path angle 
and excess (above circular) velocity, are introduced. Ifcen from the 
resulting state at the beginning of lunar parking orbit, the heights 
at pericynthion and apocynthion, y^ and Y 20 > are obtained. Now the 
state Is propagated through the time xq, resulting in a state supposedly 
directly over the landing site. Thus x 2 i and y 22 , the latitude and 
longitude of that point, are determined. The time of staying on the moon 
is used to propagate this last state to the state of the CSM at the 
time of departure from the moon’s surface. Uiis state is used to derive 
the angle y 2 ^ by which the LEM is out of the lunar orbit plane. 

From the last state, we propagate to the state at the end of lunar 
parking orbit, by introducing enough time to exhaust the total time in 
parking orbit Xy. This state is then operated on by the control para- 
meters for transearth injection and the mass y 2 ^. 

The state after transearth injection is then propagated back 
toward the earth, obtaining the state at perigee. The transearth flight 
time y 2 ^ is the time to perigee. In approximate calculations, if the 
orbit is elliptical, but very close to the center of the earth, this 
time could be negative. To create smoother convergence, an orbital 
period is added to the time in this case. At perigee, we have the return 
inclination y 2 y The next quantities are dependent on the state at 
reentry. This state is defined to have a certain flight-path angle, 
which is a function of the energy of the return orbit. From perigee 
we propagate the state back to this flight-path angle. At reentry, then, 
the height, y 2 &> and the velocity, y 29, are available. 

The calculation of the rest of the constraints is always approxi- 
mate, and is based on the assumption that the motion of the spacecraft 
during reentry may be approximately represented by a circular orbit. 
Finding the landing site may be done in either of two ways, at the ' 
option of the user. First, he may specify a fixed reentry range. The 
circular orbit is calculated as to traverse this ra^e. Thus we have 
a state at landing, frcm which azimuth at landing y 2 n, total mission 
time y^, and latitude at landing y,p are easy to derive. Alternatively, 
the user may specify the longitude or the earth landing site, in which 
case an iteration is performed to choose the time elapsed so that the 
spacecraft and the landing site have the same longitude. (For partial- 
derivative computations, this time may be adjusted by one revolution, so 
as to remove discontinuities.) From the time of landing, we get azimuth 
at landing, y2 6, reentry range, yjo, total mission time, y 3, and latitude 
of landing, y$ 2 m 
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If it is desired to pick up an initial state vector, rather than 
starting at launch, this can easily be achieved by setting a switch 
to describe to the program how far into the mission the state is. If 
enough stages of the trajectory have been computed to derive all the 
dependent variables, then the remainder is by-passed. Thus, in all 
instances, only the desired parts of the trajectory are computed, so 
as to save computer time. 

PARAMETER CORRECTION SCHEME 

Generally speaking, the first guesses applied to the trajectory 
computer will not yield values of the dependent variables that satisfy 
the constraints. This program makes use of an iteration scheme to 
correct the independent variables until the constraints are, indeed, 
satisfied, 3he scheme is described in detail elsewhere (see ref. l). 

OUTPUT 

The output section provides the value of all the converged input 
variable and the values of all the dependent variables. In addition, 
for each of the states appropriate to the converged trajectory, about 
a hundred parameters are displayed. These include coordinates of the 
spacecraft and the attracting bodies relative to several reference 
system, the polar angles corresponding to these coordinates, orbital 
elements, the orbital parameters, etc. 
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APPENDIX 

AN ENCKE METHOD ADAPTED TO MISSION ANALYSIS 
1. The Standard Encke Method 

In the standard formulation of the trajectory problem, we are given 

x*= djc = _x_ + F, (l) 

dt ^ lx? 

X = (t o ) = x o x (t o ) = x o 

where x is the vector from the central body to the spacecraft, p is the 
attraction coefficient due to the central body, F is the sum of the 
other forces acting on the spacecraft, and x , £ are the initial 
position and velocity vectors. According to°the standard Encke method, 
we introduce another differential equation 


pP 

ip? 


The solution p = p (t) of this equation, with p(t Q ) = x Q and 
p(t Q ) = x 0 can be found with extreme accuracy in closed form. 


Thus we set 


x = P + * (3) 

so that, since £ = x -P, we have to solve the differential equation 



with the initial conditions 


£ (t ) = x (t ) - p(t ) = x - x = 0 
O O o 0 o 


and £ (t Q ) = x (t Q ) - P (t Q ) = x Q - x Q = 0 . 


In the region where (l) is difficult to solve, that is, near 
|x| =0, equation (h) is much easier to integrate, so that, in general, 
the same accuracy can be obtained with less computing time. 
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2. Change of Independent Variable 

In the construction of the closed-form solution for (2), a para- 
meter 3 arises, related to t by the equation 

.t dt 

L T (5) 

r o 

In terms of 3, Kepler's equation takes the form 

t - to + - L ^ L (6) 

where f is a transcendental function of 3, and is obtained by summing 
several power series. 

If t is taken as the independent variable, eauation. (6) has to be 
solved for 3 by an iterative method, requiring numerous time consuming 
evaluation of the function f for each integration step. Using 3 as the 
independent variable, however, only requires a single evaluation. 

It remains, of course, to see what becomes of equations (2) and (4) 
if 3 is is the independent variable. We have from (5) that 


dt p 

d 3 "A 


at any point along the solution of (2)-,»A prime (') denotes differen- 
tiation with respect to 3. Thus p = p' -5— and p' = P 75— at any 
point along the solution of (2). Thus the initial conditions become 

p(3o) = x o a nd P'(Po) = X ° — 

>f4 


when 3 0 = (3(t 0 ) = 0. 


Now the solution of (2), p and p', can be written in closed form for 
any 3. As auxiliary quantities in this solution we have |p| and 


D 


P . P 1 
I PI 


They are computed as functions of 3 before p and p' known, that 
is, with accuracy at least as good as that of p and p*. Not only are 
they needed and easy to compute, but also they have the interesting 
property that 


dt _ _p 

dp JV 


, as we saw above. 


(7) 
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and dft D 

d$2 s 7 JT 

Thus equation (2) Is solved more economically in terms of 3 than In 
terms of t. 


Now we turn to equation (4) . To treat it, we want to express l " in 
terms of £ . From (7) we have that 

V = f dt = ( JLP|_ 

dp /u 

Differentiating with respect to 0, 



Thus (8) is the equation to he integrated numerical. 1 ”, instead of (4). 

The coefficients and can he calculated with much more 

P JPI 

accuracy than the factors involving $, since they depend only on the 
two-hody solution. 


For analysis of error propagation, we write (8) as 

r = — f(p+*) 4^3 - p] + F + r — * 

IP-I L IP* ^ J P IPI 


The mechanics of the procedure, then, are easy to enumerate. The 

initial conditions are x and x • Let 

o o 

pog -x o 


p* (t D ) 
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I pj d * 

Using these initial conditions, evaluate t, , ^ , P, P for each 

value of 3 to be considered. ^ 

Let = 0. Using these initial conditions, integrate equation 

(8) to get 4 (0) and (ft). Note that the first two terms on the 
right-hand side of equation (7) are functions of x and possibly x f . 
These are obtained by 

X (0) = 0(0) + £(0) 

*’(0) = p'(0) + £' (0). 

If, at any point $ is required, it can be found from 
x (t(0)) = x' (I)) 

Depending on the rectification control logic, there will be places 
where the solution to equation (2) must be started over. 

At this point, 0, I, and £' are reset to zero, while the value 
t, x, x became the new t Q , sc Q , x 0 . In particular, then, 

V = y— ie* — 

|X| 



